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The evolution of electron conductance in the presence of inelastic effects is studied as an atomic 
gold contact is formed evolving from a low-conductance regime (tunneling) to a high- conductance 
regime (contact). In order to characterize each regime, we perform density functional theory (DFT) 
calculations to study the geometric and electronic structures, together with the strength of the 
atomic bonds and the associated vibrational frequencies. The conductance is calculated by first 
evaluating the transmission of electrons through the system, and secondly by calculating the con- 
ductance change due to the excitation of vibrations. As found in previous studies [Paulsson et al. , 
Phys. Rev. B. 72, 201101(R) (2005)] the change in conductance due to inelastic effects permits to 
characterize the crossover from tunneling to contact. The most notorious effect being the crossover 
from an increase in conductance in the tunneling regime to a decrease in conductance in the contact 
regime when the bias voltage matches a vibrational threshold. Our DFT-based calculations actually 
show that the effect of vibrational modes in electron conductance is rather complex, in particular 
when modes localized in the contact region are permitted to extend into the electrodes. As an 
example, we find that certain modes can give rise to decreases in conductance when in the tunneling 
regime, opposite to the above mentioned result. Whereas details in the inelastic spectrum depend on 
the size of the vibrational region, we show that the overall change in conductance is quantitatively 
well approximated by the simplest calculation where only the apex atoms are allowed to vibrate. 
Our study is completed by the application of a simplified model where the relevant parameters are 
obtained from the above DFT-based calculations. 



PACS numbers: 72.10.-d, 73.40.Jn, 63.22.-|-m 

I. INTRODUCTION 

Recent experimental advances have permitted to 
probe electron transport processes at the atomic 
scaleji Junctions can be formed that support current 
flow through atom-sized constrictions or even single 
molecules. Atomic vibrations become detectable and 
very dependable on the environmental temperature. Ac- 
cording to the distance between electrodes, the conduc- 
tance can vary several orders of magnitude when the ap- 
plied voltages are small, typically below the cV scale. 
This behavior is due to the exponential dependence of 
current with distance when the conductance is due to an 
electron tunneling process. However, at short electrode 
distances, the current levels off and saturates: the con- 
tact regime is reached. The conductance is maximum in 
this case and a high-conductance regime is attained. The 
physics in these two regimes can be very different. 

The low-conductance regime has been thoroughly 
studied with the scanning tunneling microscope (STM). 
The initial inelastic effects were realized by showing the 
increase in conductance on an acetylene molecule when 
the bias voltage matched the C-H stretch mode<^ The 
proof that the mode was indeed excited was the iso- 
topical effect that the changes of conductance showed 
when replacing C2H2 by C2D2. This finding paved the 
way to vibrational spectroscopy with sub- Angstrom spa- 
tial resolution, permitting the identification of chemical 
components of matter on the atomic scalci^i^ The first 



experimental evidence of mode excitation in the high- 
conductance regime was achieved in monatomic gold 
wiresi^ The conductance of the wires showed clear re- 
ductions at thresholds that were proven to originate in 
the backscattering of electrons from some selected vi- 
brations of the wiresi^i^ Similarly, experiments with the 
break junction geometry have also revealed signatures in 
the conductance related to several vibrational modes of 
a single H2 molecule trapped between the electrodes.— 

The emerging picture is that in the tunneling or low- 
conductance regime, the excitation of vibrations leads 
to increases in conductance at the corresponding volt- 
age thresholds, while in the contact or high-conductance 
regime, the effect of vibrations is to reduce the con- 
ductance. Theoretical studies in the weak electron- 
vibration coupling regime have shown that the lowest 
order expansion^ is capable of correlating this behavior 
with a single parameter: the eigenchannel transmission 
probability ri^ i^^d-*^ In the simplified case of a single elec- 
tronic level connected with two electrodes under symmet- 
rical conditions, the inelastic effects (of a vibrationally 
mediated on-site modulation) go from increases in the 
conductance for r < 1/2 to decreases for r > 1/2. In 
this way, the behavior of the inelastic conductance would 
define the crossover from tunneling to contact. There is 
experimental evidence showing that this picture is indeed 
more complex. The excitation of the 0-0 stretch mode 
of the chemisorbed O2 molecule on Ag(llO)^ leads to a 
decrease of the tunneling current (instead of an increase) 



in opposition with most cases in the low-conductance 
regime ji^ii^ 

The aim of the present work is to analyze the con- 
tinuous evolution from tunneling to contact in a model 
system constituted by a junction of gold atoms, which 
provides an almost perfect realization of a single trans- 
mission channel system. The definition of when a given 
atomic system correspond to one of the two cases an- 
alyzed above is already problematic, hence we address 
this issue by investigating the behavior of different prop- 
erties of the junction with the interatomic distance. Ini- 
tially, we are interested in studying the crossover from 
tunneling to contact by evaluating the total energy, the 
strain, and the modification of vibrational modes as the 
electrode distance decreases. This allows us to find a 
range of distances where the junction behaves as either 
two independent systems or a strongly coupled one. The 
second part of this work evaluates the effect of the in- 
teratomic distance in electron transmission; this allows 
us to study the elastic conductance within Landauer's 
formalism. The correlation of the transmission against 
the interatomic relaxation permits a clear identification 
of both regimes as well as the transition region. Finally, 
the inelastic properties of the conductance are studied 
in the different regimes. The inelastic signals are inter- 
preted in a simplified model that captures the calculated 
behavior and illustrates the fundamental concepts. 

The continuous transition from tunneling to contact 
is experimentally challenging, since most metallic point 
contacts (including Au) usually exhibits a sudden jump 
in the conductance when the surfaces are brought into 
contactji^ On the other hand, experiments with a low- 
temperature STM on Cu(lll) and Ag(lll) surfaces have 
shown that both sharp jumps as well as smooth variations 
can be obtained in the crossover from tunneling to con- 
tact: when the tip is approached over a clean surface one 
observes a jump in conductance (related to the transfer 
of the tip-atom to the surface) , whereas over an isolated 
metallic adatom the evolution is smooth and reversible 
To our knowledge there is no measurement of the evolu- 
tion of the inelastic signals in the formation of a metal- 
lic point contact, likely owing the relative weak effect 
(conductance changes expected to be less than 1%) and 
difficulties with noise. Despite this we envision that our 
idealized model system is not unrealistic, and can provide 
a useful framework for investigating the complicated in- 
terplay between chemical bonding, electron conduction, 
atomic vibrations, etc. Our first-principles treatment fur- 
ther addresses all of these issues in a unified way to pro- 
vide quantitative predictions. 



II. THEORY 

The present work can be divided by the different meth- 
ods that we have used. In order to study the struc- 
tural properties of the atomic junction the standard DFT 
SiESTAfii method is used. The elastic conductance is 
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FIG. 1: (Color online) Generic setup for the calculation of 
structural properties of the atomic gold junction. The peri- 
odic supercell consists of a 4 x 4 representation of two Au(lOO) 
surfaces sandwiching two pyramids pointing towards each 
other. The characteristic electrode separation L is measured 
between the second-topmost surface layers, since the surface 
layer itself is relaxed and hence deviates on the decimals from 
the bulk values. The interatomic distance between the apex 
atoms is denoted d. 



evaluated from the transmission function of the atomic 
junction calculated with TRANSiESTAj^^ and the inelas- 
tic contribution is performed using the method presented 
in Ref. 

The system representing the atomic junction is de- 
picted in Fig. [TJ We consider a periodic supercell with a 
4x4 representation of two Au(lOO) surfaces sandwiching 
two pyramids pointing towards each other. The charac- 
teristic electrode separation L will be measured between 
the second-topmost surface layers, since the surface layer 
itself is relaxed and hence deviates on the decimals from 
the bulk values. The corresponding calculations with 
the Siesta method are carried out using a single zeta 
plus polarization (SZP) basis with a confining energy 
of 0.01 Ry (corresponding to the 5c? and 6(s,p) states 
of a free Au atom), the generalized gradient approxi- 
mation (GGA) for the exchange-correlation functional, 
a cutoff energy of 200 Ry for the real-space grid integra- 
tions, and the F-point approximation for the sampling 
of the three-dimensional Brillouin zone. The interaction 
between the valence electrons and the ionic cores are de- 
scribed by a standard norm-conserving TrouUier-Martins 
pseudo-potential generated from a relativistic atomic cal- 
culation. 

The calculations of the vibrations are performed by 
diagonalization of the dynamical matrix extracted from 
finite differences (with corrections for the egg-box effect, 
i.e., the movement of basis orbitals — following the dis- 
placed atoms — with respect to the real space integration 
grid)!^^ As the active atoms we consider initially — for 
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pedagogical purposes — just the two apex atoms and com- 
pare afterwards the results when the vibrational region 
is enlarged. 

The transport calculation naturally considers infinite 
electrodes by including the DFT self-energy calculated 
for infinite atomistic leads in the conduction equations 
Since we arc here interested in the low-bias regime (of the 
order of the vibrational frequencies) it suffices to calcu- 
late the electronic structure in equilibrium in order to de- 
scribe the elastic transport properties. While the trans- 
mission generally involves a sampling over k-points we 
here approximate it with its F-point value; this has pre- 
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viously been shown to be a reasonable approximation for 
supercells of similar dimensions in the case of atomic gold 



wires 
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Finally, the inelastic transport calculations are per- 
formed using the nonequilibrium Green's function 
(NEGF) formalism combined with the electrode cou- 
plings Tl^h extracted from the TranSiesta calcula- 
tions and the electron-vibration couplings (corre- 
sponding to modes A with energies huj\) from the finite- 
difference method.— According to the lowest order ex- 
pansion (LOE)^iiS the inelastic current reads 
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where Gq = 2e^/h is the conductance quantum, V the external bias voltage, {nx) the occupation of mode A, nF(£) 
the Fermi function, Ti. the Hilbert transform, and (3 = l/ksT the inverse temperature. The retarded Green's function 
G, the spectral function A = i(G — G^), as well as the electrode couplings Tl,r are all evaluated at the Fermi energy 
in the LOE scheme. For convenience we have also defined the quantities A^^/j = GFl^/jG^ such that A = A^ -I- A^,. 
The sums in Eq. ^ runs over all modes A in the vibrational region. For a symmetric system (such as the present one 
for the atomic junction) it can be shown that the asymmetric terms in the current expression vanishes. Furthermore, 
at low temperatures {(3 — > oo) and in the externally damped limit {{nx) ~ 0) the inelastic conductance change from 
each mode A (beyond the threshold voltage eV > hujx) is given by 
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From this expression we note that 5Gx can be either pos- 
itive or negative, depending on sign of the trace. 



III. STRUCTURAL AND VIBRATIONAL 
PROPERTIES OF THE ATOMIC JUNCTION 

As the electrode separation is decreased, we relax in 
each step the apex atoms, the base atoms of the pyra- 
mids, and the first-layer atoms until residual forces are 
smaller than 0.02 eV/A. This allows us to obtain the evo- 
lution of the (Kohn-Sham) total energy E of the system 
as a function of the electrode distance, see Fig. (5] We 
find that the energy is reduced (of the order 1 eV) by the 
attractive interaction between the apex atoms, due to the 
formation of a covalcnt bond at short distances, Fig.[51Ja). 
The slope of the energy presents a rapid change for dis- 



tances shorter than L = 16.0 A. This is more clearly seen 
in the lower part of Fig. [2] where the strain — or force on 
the unit cell — is represented. This force is evaluated as 
the numerical derivative of the total energy with respect 
to electrode separation. Here, the onset of chemical in- 
teractions is clearly seen around L = 16.0 A, Fig. [21(b); 
where the force experiences a significant increase reach- 
ing a maximum at L = 15.6 A. 

The evolution of the interaction between the apex 
atoms with distance is also revealed in the study of the 
vibrational modes. This is presented in Fig.[3l where the 
blue connected data points correspond to the 6 modes 
where only the apex atoms vibrate, and the red crosses 
to the 30 modes where also the pyramid bases vibrate. 
These modes follow different behavior with the electrode 
separation. 

In the following we analyze the simplest case of just 
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Electrode separation L [A] 

FIG. 2: (Color online) Total energy differences and the nu- 
merical derivatives as a function of the electrode separation. 
The lower part of the figure describes the strain on the unit 
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FIG. 3: (Color online) Vibrational frequencies versus elec- 
trode displacement. The connected data series refer to the 
situation where only the two apex atoms are vibrating (result- 
ing in the six vibrational modes indicated in the plot); circles 
symbolize the two longitudinal modes (CM and ABL) and 
diamonds the four (pair-wise degenerate) transversal modes. 
The red stars are the corresponding vibrational frequencies 
when also the pyramid bases are considered active. The three 
regimes are clearly identifiable: (a) concerted apex vibrations, 

(b) crossover where the stretch modes become degenerate, and 

(c) independent apex vibrations. 



the two apex atoms. Generally, two longitudinal stretch 
modes (represented with connected circles in Fig. [3|) line 
up the highest in energy. For an electrode distance larger 
than L = 16.5 A these correspond to the isolated (i.e., 
decoupled and hence degenerate) stretch modes of each 
apex atom. Fig. [3l^c). As the electrodes are approached, 
the attractive apex-apex interaction leads to a slight dis- 
placement of the apex atoms away from the base of the 
pyramids. The consequence is a small weakening of the 



apex-atom coupling to the base which results in decreas- 
ing frequencies, i.e., to softening of the modes. Another 
consequence of the increasing interaction is the splitting 
of the degenerate modes into a symmetric (out-of-phase) 
and an antisymmetric (in-phase) mode. We will refer to 
these as the alternating bond length (ABL) mode^ and 
the center of mass (CM) mode, respectively. When the 
electrode separation reaches the region between L = 15.8 
A and L = 16.0 A the frequencies drop significantly. 
Fig. EKb). This points again at the chemical interaction 
crossover that we presented in the previous paragraph: 
now the interaction between the apex atoms becomes 
comparable with the interaction with the electrodes and 
hence weakens the stretch modes initially set by the in- 
teraction between the apex atom with the base of the 
pyramid. As the apex-apex interaction grows larger, the 
modes start to increase in frequency and further show an 
significant split, Fig.[3]Ja). 

The behavior of the two stretch modes of Fig. [3] is 
easily understood with a simple one-dimensional elastic 
model of two masses, each coupled to infinite-mass sys- 
tem with a spring constant ki, and interconnected by 
another spring constant k2. The frequencies of the two 
stretch modes are then wcm = \/ki/m (in-phase) and 
i-^ABL = \J (ki + 2fc2)/m (out-of-phase) , where m is the 
mass of each atom. Note that in the tunneling regime 
the apex-apex interaction is attractive, cf. Fig. [21 which 
would correspond to a negative value of ^2. When the 
bond has been formed, k2 can be represented as classical 
(positive) spring constant. This model essentially cap- 
tures the evolution of the stretch modes. In particular, 
the sign change of fc2 a-t the chemical instability explains 
the mode crossing between L = 15.8 A and L = 16.0 A, 
Fig. ^h), and why the CM mode has a higher frequency 
than the ABL mode in the tunneling case, and vice versa 
in the contact case. 

The analysis of the modes with electrode distance thus 
permits us to recover the same range of distances with the 
chemical crossover as in the preceding section concerning 
the total energy and strain. This identification is also 
possible from the more realistic calculation that includes 
the vibration of the base atoms. 



IV. ELASTIC CONDUCTANCE 

In this study we are interested in the low-bias regime. 
Hence the elastic conductance is determined via Lan- 
dauer's formula by the transmission at the Fermi energy 
ep. As expected for the gold contact, wc find that the to- 
tal transmission is essentially due to a single eigenchannel 
(for the geometries considered here the contribution from 
the secondary channel is at least three orders of magni- 
tude smaller). Figure 3] plots the transmission r and the 
apex-apex distance d as a function of electrode separa- 
tion L. In the tunneling regime the transmission is char- 
acterized by an exponential decay with separation. It is 
instructive to compare this with the transmission proba- 
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Electrode separation L [A] 

FIG. 4: (Color online) Transmission r (blue disks) and apex- 
apex distance d (red crosses) versus electrode separation L. 
In the tunneling regime the transmission decays exponentially 
with separation as indicated with the dashed line (correspond- 
ing to a tunneling parameter A = 1.54 A ). The point at (a) 
corresponds well with the contact region of transmission one 
and closest apex separation, (b) is near half transmission and 
the instability in apex separation, (c) is finally the tunneling 
regime, where the apex atoms are independent. 



bility T cx exp(— 2Ai5) for a rectangular one-dimcnsional 
barrier, where A = is a characteristic tunnel- 

ing length, $ the apparent barrier height, and D the 
barrier width (valid for AD ^1). As shown in Fig. 21 an 
exponential fit to the calculated tunneling data leads to 
a tunneling parameter A = 1.54 A which would corre- 
spond to an apparent barrier height of the order $ « 9.1 
eV. Compared with measurements of the work function 
on perfectly flat Au surfaces (5.31-5.47 eV)^ this value 
is certainly high. On the other hand $ is not very well de- 
termined from an exponential fit to data spanning only 
one decade. Another contribution to a relatively large 
barrier height could be geometric effects from the pyra- 
midal shapes. 

The deviation from the exponential tunneling behavior 
(visible around L — 16.0 A) is a clear indication of the 
crossover to contact. The contact regime is characterized 
by a constant transmission equal to unity since an atomic 
gold junction has effectively only one conduction chan- 
nel. The value t = 1/2 to define the crossover between 
contact and tunneling is somewhat arbitrary, however a 
detailed comparison with the previous section justifies 
this definition. Indeed, Fig. [5] also shows the behavior of 
the apex-atom distance d with electrode separation, per- 
mitting us to make contact with the chemical crossover 
defined previously. Between L = 15.8 A and 16.0 A we 
find that the apex-apex distance has changed by almost 
0.7 A. This shows that at these electrode distances, there 
is an instability that drive the formation of a covalent 
bond between apex atoms. This agrees with the conclu- 
sion from both total energy, strain and frequency calcu- 
lations that the crossover takes place between 15.8 A and 
16.0 A. From Fig. 2] we identify a transmission 1/2 asso- 



ciated with i = 15.9A(d = 3.7A) hence permitting us 
to identify the crossover from tunneling to contact with 
the chemical crossover. 



V. INELASTIC CONDUCTANCE 

The behavior of the inelastic contributions to conduc- 
tance is very different in the two studied regimes. In 
the tunneling regime the opening of the inelastic channel 
enhances the conductance of the system, while the cre- 
ation of a vibrational excitation in a high conductance 
regime is a source of backscattering that decreases the 
conductance. Figure [5] shows the calculated change in 
conductance (second derivative of the current with re- 
spect to bias voltage cPl/dV^) for the contact, crossover 
and tunneling regions. These three typical cases — labeled 
(a), (b), and (c), respectively — are indicated in the pre- 
vious Figs. [2H4] for easy reference. We investigate how 
the inelastic conductance change depends on how many 
atoms in the junction that are considered active: in Fig. [5] 
the thick blue line is the spectrum corresponding to only 
the two apex atoms vibrating, the dashed red curve to 
the 10 pyramid atoms vibrating, and the dotted black 
curve to the pyramids and first-layer atoms vibrating (42 
atoms). In this way we follow the convergence of the cal- 
culations as the vibrational region is gradually enlarged. 
The essential data from these calculations are summa- 
rized in Tab. [J 

From the simplest case when only the two apex atoms 
are vibrating, we arrive at the conclusion that only the 
two longitudinal stretch modes contribute to the change 
in conductance, leading to the qualitative known result 
of increase of the conductance in tunneling regime and 
decrease in contact. The crossover case Fig.[5l[b) presents 
a combination of an increase in conductance from the 
ABL mode and a decrease from the CM mode. 

This behavior is a signature of the different processes of 
conduction. In the tunneling case, the tunneling process 
is determined by the more slowly-decaying components 
of the electron wave function of the surface. Because 
of the exponential tunneling probability dependence on 
distance a mode that modulates the tunneling gap is ex- 
pected to contribute positively to the conductancci^ In- 
deed this is the case for the ABL mode. Furthermore, 
the CM mode that correspond to a fixed apex-apex dis- 
tance cannot contribute positively, neither the transverse 
modes because none of them decrease the apex-apex dis- 
tance from the equilibrium position during a vibration 
period. Instead, the CM mode is found to contribute 
negatively to the conductance, cf. Tab. U A simplified 
model presented in the next section will explain this ob- 
servation. 

In the contact case, the electronic structure responsible 
for the conduction process is largely concentrated upon 
the apex atom, hence the transport is being modified by 
the motion of basically only these atoms. Indeed both 
the ABL and CM modes lead to drops in the conductance 
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FIG. 5: (Color online) Second derivative of the current versus bias voltage for three characteristic situations (a) contact, (b) 
crossover, and (c) tunneling. In each situation we consider different active vibrational regions: the two apex atoms only (thick 
blue line), the 10 pyramid atoms (thick dashed red curve), and both pyramids and first-layer atoms (dotted thin black curve). 
The signal broadening is due to temperature {T = 4.2 K). 



as is evident from Fig. [5]Ja) and Tab. [H The transverse 
modes give essentially no signal; this is similar to the 
findings for atomic gold wires where the transverse modes 
cannot couple because of symmetryi^ii^ 

Figure [5] shows how the inelastic spectrum is modified 
if wc increase the vibrational region by allowing more 
atoms to vibrate. In the tunneling and contact cases we 
sec that the single main peak splits up into a number of 
peaks, indicating that the apex vibrations are actually 
coupled with the vibrations in the bulk. For the contact 
case the broadening of the signals is expected to be di- 
rectly influenced by the phonon density of states of the 
bulk. As was first shown by Yanson^S^ the spectroscopy 
of microcontacts at low temperatures — a technique nowa- 
days referred to as point contact spectroscopy — reveals 
a signal in (Pl/dV^ which is a direct measurement of 
the Eliashberg function a^F, i.e., roughly speaking the 
product of the squared electron-phonon coupling matrix 
element a and the phonon density of states F, averaged 
over the Fermi sphere.— In the case of microcontacts the 
measured signal is predominantly due to the transverse 
modes. This is in contrast to our case of the atomic point 
contact, where wc only find signals from the longitudi- 
nal modes. However, from Fig. [S{a) we notice a signal 
broadening by increasing the vibrational region, pointing 
towards the vibrational coupling to the bulk modes. 

In the crossover region between tunneling and contact, 
Fig. \^h) shows a dramatic change depending on the size 
of the vibrational region. Different modes give positive or 
negative contributions in the conductance, but in such a 
way that they lead to an overall absence of (or relatively 



small) variation in the conductance, cf. Tab. HI 

Comparing the total change in conductance 
AG = G{V > hwx) - G{V = 0) induced by all modes 
(for the tunneling, crossover, and contact situations), 
wc find that the calculations with different vibrational 
regions give almost the same results. As found in 
Tab. [H we thus conclude that to a first approximation 
we can describe AG « SG abl + SGqm, i-e., the overall 
conductance change can be estimated with the minimal 
vibrational region (the two apex atoms). This simple 
approach does however not accurately describe details 
of the (Pl/dV"^ spectrum. 

VI. DISCUSSION 

The effect of the tunneling to contact crossover has im- 
portant implications in the inelastic conductance since in 
the first case the inelastic effects trend to increase and 
in the second case to diminish the electron conduction. 
From the results of the previous section, we have seen 
that this crossover roughly takes place at the same range 
of distances as in the elastic conductance case. By look- 
ing at the transmission in the elastic conductance case, 
we conclude that when the transmission is t = 1/2 we 
should also be near the crossover between tunneling to 
contact in the inelastic case one. This finding is simi- 
lar to the crossover found for the single-state impurity 
model analyzed in Ref. 0]. However, in the present case, 
the system is not obviously modeled with a single-state 
impurity. Instead we can easily reproduce the same kind 
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"Only apex atoms vibrating, device includes first-layer atoms 
'Apex and base atoms vibrating, device includes first-layer atoms 
'^Pyramids and first-layer atoms vibrating, device includes first- 
and second-layer atoms 

TABLE I: Characteristic data for the six structures representing the evolution of the junction from tunneling to contact regimes. 
The columns describe the electrode separation L, apex-apex distance d, transmission t, vibrational energies fiujx and mode- 
specific conductance changes SG\ (for the ABL and CM modes), and the total conductance change AG/G from all modes 
(calculated for three different sizes of the vibrational region). 



of analysis for a slightly more sophisticated model, where 
two impurities are connected to reservoirs and interacts 
via a hopping term between them. Under symmetric con- 
ditions this system is described by 



H 



£0 t 




7 













t £o 
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(5) 



where the Hamiltonian H includes on-site energies and 
a hopping matrix element t. The level broadening func- 
tions Ta describes the coupling of the sites to the con- 
tacts a = L, R with strength 7 (which in the wide band 
approximation are considered energy independent). The 
corresponding retarded Green's function is 



G 



H 

2 



iiTi 



(2Ae + i7)2 - 4t2 



2Ae + Z7 2t 
2t 2Ae + 



(6) 



where in our case Ae = ep — eo ^ 7 holds since the level 
positions would be close to the Fermi energy ep (the on- 
resonance case). The transmission becomes 



16^7 



2^,2 



(4t2 + 72)^ 



0{As% 



(7) 



where perfect transmission r = 1 is obtained for t ~ 7/2. 

Inspired by our electron-phonon coupling matrices ob- 
tained from the full DFT calculations, we assign the fol- 
lowing forms to the longitudinal ABL and CM mode cou- 
plings 



Ma: 



BL 



mi m2 
m2 mi 



M 



CM 



TO3 









(8) 



The ABL mode is symmetric and generally described by 
two coupling strengths: mi represents an on-site modifi- 
cation via a change in the electrode coupling, whereas 1112 



is a modulation of the hopping between the apexes. Cor- 
respondingly, the CM mode which is asymmetric bears 
an asymmetric on-site modulation and no hopping 
modulation (fixed apex-apex distance). With these ex- 
pressions we can simply evaluate Eq. ^ to find the fol- 
lowing inelastic conductance changes 

{l6t^ - 2472t2 + 74) 



SG 



ABL 



Got 



(4*2+72)^ 
2 16miTO2Ae 



m2 



SG. 



CM 



Got 



t2 

16m|72 

" (4i2 + 72)2 



4^3 



0{Ae'), 



OiAs')], (9) 
(10) 



We first discuss the conclusions to be drawn about the 
ABL mode. Notice that i5Gabl is only weakly dependent 
on the the on-site coupling element mi and vanishes on 
resonance (Ae = 0). In the tunneling limit {t 0) we 
find that 



lim 

t^o 



<5Gabl 



16m2 



OiAe'), 



(11) 



Go 7^ 

i.e., the ABL mode gives a positive contribution to 
the conductance proportional to the square of coupling 
strength TO2. In the contact limit (t ^ 1) we find 

SGabl 4m2 16mim2Ae 



0{Ae'), (12) 



lim 

*i Go T T 

i.e., the ABL mode gives here a negative contribution. 
The exact crossover between an increase and a decrease 
is determined by solving SG abl — 0, which indeed yields 
T = 1/2 as is the case for the single-site casci^ 

Next, we see from Eq. (fTO| that the conductance 
change i5Gcm from the CM mode is always negative (i.e., 
the CM mode backscatters even in the tunneling regime) . 

These results thus permit to rationalize the crossover 
from tunneling to contact for the inelastic conductance — 
as found numerically in Sec. |V] — as taking place around 
a transmission of t 1/2. 
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VII. SUMMARY AND CONCLUSIONS 

The evolution of the inelastic signals from the tunnel- 
ing to contact regimes has been studied through DFT 
calculations. We have compared our results with the 
crossover between bonding and rupture of the atomic 
junction by studying the geometric and electronic struc- 
tures of the junction, together with the strength of the 
atomic bonds and the associated vibrational frequencies. 
This permitted us to find a typical transition distance 
between electrodes where a small change leads to a large 
readjustment of the apex-apex atom distance, as well as 
a change of the strength of interactions as revealed by 
the total energy, the strain and the frequencies of the 
system's modes. 

The conductance has been calculated by first evaluat- 
ing the transmission of electrons through the system, and 
second by calculating the conductance change due to the 
excitation of vibrations. As found in previous studies^ 
the change in conductance due to inelastic effects permits 
to characterize the crossover from tunneling to contact. 
The most notorious effect being a decrease of conduc- 
tance in the contact regime to an increase in the tunneling 
one when the bias voltage exceeds the vibrational thresh- 
olds. Our DFT-based calculations show that the effect of 
vibrational modes in the cPl/dV'^ spectra is rather com- 
plex, in particular when modes localized in the contact 
region are permitted to extend into the electrodes. As 
an example, we find that certain modes can give rise to 
decreases in conductance when in the tunneling regime, 



opposite to the above mentioned result. Whereas details 
in the inelastic spectrum depends sensitively on the size 
of the vibrational region, we find that the magnitude of 
the overall change in conductance can actually be reason- 
ably described with just the minimal case where only the 
apex atoms vibrate. This means that while the modes 
are rather delocalized the region of inelastic scattering is 
localized around the apex atoms. 

By comparing our results with a simplified model, we 
conclude that in this single eigenchannel problem the tun- 
neling to contact crossover takes place exactly at t = 1/2, 
in agreement with the findings for the elastic conduction 
process and the chemical crossover. Hence, we can trace 
back the origin of the conduction process, both in the 
presence and absence of vibrational excitation, to the 
same type of underlying electron structure that deter- 
mine the electrode's chemical interaction and the electron 
conductancei^ 
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